Steady states of a x 3 parametric oscillator with coupled polarisations 
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Polarisation effects in the microcavity parametric oscillator are studied using a simple model 
in which two x 3 optical parametric oscillators are coupled together. It is found that there are, 
in general, a number of steady states of the model under continuous pumping. There are both 
continuous and discontinuous thresholds, at which new steady-states appear as the driving intensity 
is increased: at the continuous thresholds, the new state has zero output intensity, whereas at 
the discontinuous threshold it has a finite output intensity. The discontinuous thresholds have no 
analog in the uncoupled device. The coupling also generates rotations of the linear polarisation of 
the output compared with the pump, and shifts in the output frequencies as the driving polarisation 
or intensity is varied. For large ratios of the interaction between polarisations to the interaction 
within polarisations, of the order of 5, one of the thresholds has its lowest value when the pump is 
elliptically polarised. This is consistent with recent experiments in which the maximum output was 
achieved with an elliptically polarised pump. 



PACS numbers: 42.65. Yj, 71.36.+C, 42.25. Ja, 78.20.Bh 
I. INTRODUCTION 

Semiconductor microcavities are high finesse Fabry- 
Perot structures, typically consisting of a planar semi- 
conductor cavity layer bounded by Bragg mirrors. The 
mirrors confine two-dimensional photons, which mix with 
the exciton states of quantum wells embedded in the cav- 
ity. Such mixing gives a type of two-dimensional polari- 
ton known as a "cavity polariton'—. The non-linear dy- 
namics of cavity polaritons corresponding to polariton- 
polariton scattering has been extensively studied, due to 
the possibility of observing bosonic effects such as stimu- 
lated scattering. Recent experiments have demonstrated 
new aspects to the non-linear dynamics of coherent po- 
laritons: parametric oscillation and amplification. 

Parametric oscillation and amplification in microcav- 
ities has been demonstrated using both pulsed^ and 
continuous- wave* excitation. In these experiments, a 
laser tuned near to the energy of the lower polariton mode 
is used to generate a coherent polariton field in the mi- 
crocavity. Owing to a x^ non-linearity provided by the 
exciton-exciton interaction, this pump mode is coupled 
to "signal" and "idler" modes at lower and higher ener- 
gies. The coupling corresponds to the scattering of pairs 
of pump polaritons into the signal and idler, in contrast 
to the coupling in the conventional x^ optical paramet- 
ric oscillator—, which corresponds to the fission of pump 
photons. Above a critical pump intensity, the gain due to 
the nonlincarity outweighs the damping of the signal and 
idler modes. When this occurs, the steady-state in which 
there is a single coherent field at the pump becomes un- 
stable towards a state which also has coherent fields at 
the signal and idler. In the continuously pumped exper- 
iments, this instability develops spontaneously, and the 
new steady-state is reached. In the pulsed experiments 



ofRef. H however, there is not enough time for the in- 
stability to develop spontaneously before the excitation 
pulse is over. Instead it is triggered using a second "seed" 
laser pulse injected into the signal mode. 

The theory of microcavity parametric oscillation and 
amplification was initially developed by Ciuti et al& and 
Whittakeri In the former, the pulsed measurements of 
Ref. are treated within a quantum optics formalism, 
while the latter used classical nonlinear optics to explain 
the steady-state behaviour. The dynamical equations 
which occur in both models are the same, demonstrat- 
ing that the phenomena are essentially classical effects, 
with the exception of the incoherent luminescence which 
occurs below threshold— More recent theoretical work 
by Savasta et al— includes frequency-dependent nonlin- 
earities and nonlinear absorption, arising from exciton- 
exciton correlations. 

Parametric oscillation requires a significant electro- 
magnetic response at the wavevectors and frequencies of 
the signal and idler. Thus the signal and idler must lie 
near the polariton dispersion. The signal and idler must 
also satisfy the requirements of wavevector and frequency 
conservation in the generation process. In the conven- 
tional parametric oscillator, these two requirements are 
usually met by exploiting birefringence^. Thus the po- 
larisations of the fields for which the device operates are 
prescribed. In the microcavity parametric oscillator how- 
ever, the unusual dispersion of cavity polaritons allows 
them to be met irrespective of the polarisations. This 
is achieved for pump fields near to a particular "magic" 
wavevector. In this paper we study the effects of these 
polarisation degrees-of- freedom. 

There are two recent experiments on the polarisation 
effects in microcavities that are resonantly pumped near 
to the magic wavevector, one using continuous wave 
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excitation^, and one using pulsed excitationAi. Both 
these papers argue that their observations imply the ex- 
istence of interactions between polaritons of different cir- 
cular polarisations. While some aspects of the pulsed 
dataii were recently explained using a model without 
such interactions^, they may still be necessary to un- 
derstand the steady-state experiments of Ref. 10. Fur- 
thermore, interactions between polaritons of different po- 
larisations seem to be necessary to explain the polari- 
sation dependence of four-wave mixing experiments in 
microcavities^. They could originate from interactions 
between excitons of different polarisations, which have 
been used to explain the polarisation dependence of four- 
wave mixing in quantum wells^*^. 

Because the parametric oscillator involves the dynam- 
ics of coherent polaritons, and not simply scattering, the 
consequences of an interaction between polaritons of dif- 
ferent polarisations are not obvious. In this paper, we 
investigate the effects of polarisation coupling on the 
steady-states of a simple model. 

The remainder of this paper is organised as follows. In 
section [n] we present the model. In section IIIII we ex- 
plore the steady-states of the model. The steady-state 
calculation is subdivided: in section IIII Al we calculate 
the possible values of the pump polariton fields, in sec- 
tion llll Bl we calculate the output fields for those values of 
the pump fields, and in section llll Cl we combine these re- 
sults with the "pump depletion" equations to determine 
the behaviour for a particular external drive. In section 
IIVI wc qualitatively compare our results with the con- 
tinuously pumped experiments reported in Ref. fToL and 
comment on the stability of our solutions and possible 
microscopic origins for our phcnomenological coupling. 
Finally, section Ivl summarises our conclusions. 



Here j p etc are the homogeneous line widths of the po- 
lariton states, and c p and the amplitudes of the 
photon and exciton in the polariton, i.e. the Hopficld 
coefficients. They appear because it is only the exci- 
tonic part of the polariton which interacts. f-\(t) is the 
external driving field for the pump, in a frame rotat- 
ing at the appropriate polariton frequency For the con- 
tinuously pumped situation we consider in the present 
paper, = /f exp(— iS p t), where S p is the pump de- 
tuning. The nonlinear term in the pump equation lla|) 
describes the scattering of pairs of polaritons out of the 
pump field, while the terms in l|lbllc|) provide the cor- 
responding growth in the signal and idler fields. The 
nonlinear exciton blue-shiftS^ is neglected for simplicity. 



To incorporate the polarisation degrees-of-frcedom, we 
introduce pump, signal and idler fields for the other cir- 
cular polarisation, denoting their exciton amplitudes by 
Pl, s^and i±. Without any coupling terms, their dynam- 
ics is given by the analogs of However, we now in- 
troduce a second excitonic nonlinearity which cou- 
ples the up and down spin excitons. Rather than derive 
a realistic microscopic model of exciton spin scattering, 
we treat this phcnomcnologically by choosing the simple 
form 5o I '/'T 1 2 1 ^1 1 2 ■ We assume that the coefficients k and 
g are independent of momentum and energy. While the 
former is physically justified because the wavelengths of 
the polaritons which scatter are much larger than any 
excitonic length scaled, it is more difficult to rule out the 
possibility of an energy dependent process. 



II. MODEL 

The model we analyse in this paper is a generalisa- 
tion of the scalar treatment described in Ref. 0. That 
model considers the scattering between pump, signal and 
idler polaritons of one circular polarisation. The exci- 
ton amplitudes in these fields are denoted by p^, s^and 
if. They are time dependent, so the exciton field at 
the pump wavevector k p , for example, takes the form 
Pf (t) expi(k p .r — to p t), where u> p is the lower branch po- 
lariton frequency. For simplicity, we assume that the bare 
polariton frequencies satisfy the triple resonance condi- 
tion 2uip — + We also assume that the time de- 
pendence of the amplitudes is slow compared with the 
polariton splitting, so the upper branch can be omitted 
from the model. The scattering is modelled by a term 
proportional to |</>f| 4 in the Lagrangian density, corre- 
sponding to a x' 3 ) nonlinearity. The equations governing 
the exciton amplitudes are then 

-^(|+7p) PT +2kWt = (la) 
{jt +ls ) Sl + KP1 ^ = ° (lb) 



The interaction between excitons of opposite spins in- 
troduces two new types of polariton scattering terms into 
the equations of motion. If g is constant as wc assume 
these both have the same strength, g , but for now we 
add subscripts in order to distinguish the processes in 
Eqs.[3 The first process we describe as cross-polarisation 
parametric scattering, where a pair of pump polaritons 
of opposite polarisation scatter into a signal and idler 
modes, also of opposite polarisation. These terms are 
written with a coefficient g c . The second process we de- 
scribe as a polarisation-flip, where a pump and signal, or 
idler, polariton exchange polarisations. This flip process 
is given a strength g s . A similar flip process can also 
occur between signal and idler polaritons, but this possi- 
bility is neglected here, making the assumption that the 
signal and idler amplitudes arc small compared with that 
of the pump. 



With these polarisation-coupling terms, the equations 
for the spin-up fields are 
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+ -f s J s T + ftpf i* + ffcPTPi*! + ffsPTPi s l = °> 



\x s \- \dt 
i / d 

—2 I ^ + 7i ] *t + npfs* + g c pwi s l + gspwln 



= 0. 



(2a) 
(2b) 
(2c) 



The equations obeyed by the spin-down fields are given 
by flipping the spin labels in Eqs.El 

In what follows we will not be concerned with the ab- 
solute intensities of the fields or external pumps. This 
allows us to eliminate the normal coupling k by scal- 
ing all the fields and the pumps according to — > 

Q 

~Xyf]^/~K = ^T> an d P\ ~ * P]/V^ °t c -- After this rescal- 
ing, the coupling strengths in J2J are replaced by their 
ratios to k: g c — ► g c / n = g etc. 



III. PARAMETRIC OSCILLATION 
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Pump-down intensity, L/I. 
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Pump-down intensity, Ij/I 



For particular amplitudes of the pump fields, equations 
12 bl and [2c| admit harmonic solutions with finite ampli- 
tudes for the signal and idler fields. To find these steady 
states, we set p-\(t) = pf exp (—iS p t) etc. When the de- 
tunings obey 

2S p = 5 a + 6i, (3) 

the equations for the signal and idler amplitudes become 
time- independent; defining complex rescaled detunings 
= -(6p + vy p ') = -(S p + h P )/\ x p\ 2 Gt c- they read 



FIG. 1: Left panel: allowed intensities of the pump fields 
giving steady-state operation with g — 0.2(solid lines), and 
g = (dashed lines). The sloping line and dots are for com- 
parison with Fig. [I] The sloping line is a circularity of a — 0.1 
for the pump fields, corresponding to the circularity of the 
drive used for Fig. 01 The dots mark the pump fields at the 
thresholds of Fig. [I] Right panel: corresponding signal detun- 
ing for g = 0.2. For both plots 7„ = 0.25meV,7i = l.OmeV, 
\x s \ 2 = 0.5, \Xi\ 2 = 0.97, and 5 P = 0. The intensities of the 
pump fields are given relative to that for steady-state opera- 
tion with a single polarisation, Iq. 



a s s t +Pi 2i * + gpiPii*i + mp*i s i = °> ( 4 ) 

AJa* +p* 2 s T +gp*p*si + gp*Pii* = 0. (5) 

Along with their spin-flipped counterparts, (@J and (0) 
form a set of linear homogeneous equations, parametrised 
by the pump amplitudes, for the signal and idler ampli- 
tudes: 



M 



= o. 



(0) 



The matrix of coefficients M, combined with the condi- 
tion Q, determines pump amplitudes and detunings for 
which steady state operation is possible, and the signal 
and idler fields in these steady states, as a function of 
7s, 7-tj S p , g and the Hopfield coefficients. 



A. Allowed Pump Fields and Detunings 

To determine the pump amplitudes and detunings for 
which steady-state operation is possible, we note that 
since (0 is homogeneous, solutions with a finite signal 
and idler are only possible if M has a zero eigenvalue. 
In terms of the pump polariton intensities ij = |p|| 2 , 



\Pl 



| , this occurs when 
(A S A*-/ T 2 )(A S A*-/, 2 ) 



g 2 hl 



(2/ T -A s -A*)(2/ i -A s 



a;), 



(7) 



which directly determines I-\ as a function of ij and the 
detunings. Since the determined I-\ should be real it also 
gives an equation, parametrised by 7| , among the detun- 
ings, which combines with Q to determine S s and Si as 
a function of 7j. 

In Fig^ we illustrate the pump fields and signal detun- 
ings giving steady-state operation, for a resonant pump 



with 7 S = 0.25mcV, ji — lmeV, g = 0.2, and the Hop- 
field coefficients \x s \ 2 = 0.5 and \Xi\ 2 = 0.97. We have 
estimated these damping rates and Hopfield coefficients 
to be those of the experiment reported in Ref. 0. The 
curves describing the allowed pump fields look very sim- 
ilar to those for the uncoupled device, shown as dashed 
lines, except that the degeneracy, where both polarisa- 
tions are on threshold, has been split. The signal de- 
tuning shows small deviations from the uncoupled case, 
where it would be zero with this resonant pump. These 
shifts of the signal detuning from its uncoupled value are 
due to the spin-flip processes and the imbalance in the 
damping of the signal and idler; without spin- flips or for 
7 S ' = ji' the condition for the intensities to be real is 
QAsA* = 0, as in the uncoupled device. 

We can determine whether the form shown in Figgis 
general for small g by using perturbation theory to cal- 
culate how the degeneracy is split by the coupling. We 
expand the left-hand side of Eq.[7]to first-order in the de- 
viation of the pump intensities from the degeneracy and 
the change in the detunings compared to the uncoupled 
case, and take the right-hand side to be unchanged to 
leading order. Eliminating the detunings such that the 
fields remain real, we find that the shifts in the pump in- 
tensities obey a real quartic form. Owing to its complex- 
ity, we have not studied this quartic in general. However, 
for the special case of a resonant pump it becomes 

Ml (a' - 4/ 2 « T ) (<J/ T + SI l ) 2 + a" 2 = 0, (8) 

where a = a' + ia" is the right-hand side of Eq. \J\ evalu- 
ated on the degeneracy. Eq. [8] always has the form seen 
near the degeneracy in Fig. ^ so the structure of that 
figure is general for resonant pumping and small g. 

We can also use perturbation theory away from the 
degeneracy, to study the shift in the allowed pump fields 
produced by a small g. Again for a resonant pump, the 
leading deviation in 7f from an uncoupled solution in 
which Ij = Iq is on threshold while ij is well away from 
it obeys 

»' = - 2( 7 ,V-/?) hvW-(v-7,y). 

Thus for a small value of one polarisation, turning on the 
cross-polarisation coupling increases the threshold for the 
other polarisation. There is therefore a region of I± , just 
above the uncoupled threshold, in which is multi- 

valued, although this is not visible on the scale of Fig. 
n Such behaviour is possible because in the uncoupled 
device the signal and idler fields of the bclow-threshold 
polarisation are zero, and hence there is no loss through 
these channels. With a finite coupling however, these 
fields become finite, providing another loss channel. The 
modes do not always mix in this way however; for general 
pump detunings the initial shift can be to higher or lower 
thresholds. 

In the strong-coupling limit, g — > oo, the solutions to 
Equation [7| for 7s 7^ 7- have either /j or I± = 0. The ap- 
proach to this limit is illustrated in FigEl where we plot 




FIG. 2: Allowed pump-up intensity as a function of pump- 
down intensity for g — 0.2(thin solid lines), 1.0(thin dashed 
line), 2.0(thin dotted line), 5.0(thick solid lines), 7.0(thick 
dashed lines), 10.0(thick dotted lines). The remaining pa- 
rameters are as used for figure 

the allowed pump fields for increasing values of g. The 
lower branch simply collapses towards the origin. The 
upper branch disappears, then reappears as two disjoint 
branches with asymptotes Ij = and I± — 0. These 
branches then merge, before finally collapsing into the 
origin. 

B. Signal and idler fields 

We now consider the signal and idler fields in the 
steady-states. These fields are determined, up to an over- 
all complex scale factor z, by the eigenvector of M cor- 
responding to the zero eigenvalue. 

The eigenvectors of M, unlike the eigenvalues, depend 
on the phases of the pump fields, <^>f and </>^. This depen- 
dence can be extracted by noting that M can be written 
in the form S'MS, where S is a diagonal matrix with 
entries e"^ ; e ^r , e ~^i ; . Thus the phases of the 
pump fields simply shift the arguments of the steady- 
state signal and idler fields: supposing e= (s-|-, i*, s±, i|) 
is an eigenvector of M when <f>^ = (j>± = 0, then 

d = (sTe^Mte-^s^MJe-^) is the correspond- 
ing eigenvector for finite phases. The phases of the pump 
fields have such a simple effect because there is no phase 
dependence in the form of the interaction energy we have 
chosen. A nonlinearity such as (5i</>) 4 , rather than |c/>| 4 , 
might lead to a more complicated effect of the pump 
phases. 

In Fig. |2l we plot the components of the eigenvector of 
zero eigenvalue for the pump fields shown in Fig. ^ We 
have normalised the eigenvector so that the total inten- 
sity is one and the phase of the signal up field is zero, 
and taken the pump fields to be real and positive. The 



FIG. 3: Signal fields(solid lines) and conjugate of the idler 
fields (dashed lines) corresponding to the steady-states shown 
in Fig. Q scaled such that the total intensity is one and the 
phase of the signal up component is zero. 



idler fields are always smaller than the corresponding sig- 
nal fields due to the stronger damping of the idlers. The 
crossings of the signal curves occur when the intensities 
in the two pump components are equal. In the region 
I[/Iq < 1, the polarisation with the largest pump fields 
also has the largest signal and idler fields, but that order- 
ing is reversed in the region Ijh > I. For IJIq < 1, the 
phase differences between corresponding components in 
the two polarisations lie near to zero, while for if /io > 1 
they lie near to ir. 

The two circular components of the signal field can 
be combined to form, in general, an clliptically polarised 
state. The phase differences between the two components 
of the signal that can be seen in Fig|3| correspond to ro- 
tations of the ellipse describing the signal polarisation 
compared with that describing the pump polarisation. 
Such polarisation rotations are absent in our model if 
there are no spin-flip processes. 



C. Dependence on the driving fields 

The steady-state reached in the device is selected from 
the possibilities shown in Fig. ^ by the external driving 
fields, according to the steady-state version of the pump 
equation l|2a|) 

A p p T + 2s T i T p| + g(sfif + sthH ( 9 ) 
+g(s r s* l +iii* L )p l = F h 

and its spin flipped counterpart. The first term on the 
left-hand side of Eq. describes the bare response of 
the pump field, while the remaining "pump depletion" 
terms describe the effect on the pump fields of the non- 
linear processes which generate the signal and idler. The 
pump equations 10 determine the remaining four real 



12 3 4 

Driving intensity 

FIG. 4: Steady-states, for the parameters of Fig. Q with an 
external drive of circularity a = 0.1 and varying intensity. 
The intensity of the drive is expressed relative to that of the 
lowest threshold in the uncoupled device with the same pa- 
rameters. 

unknowns: 7j_, the arguments of the pump fields, and the 
total intensity of the output fields, \z\ 2 . The pump equa- 
tions are independent of the overall phase of the zero 
eigenvector, argz, corresponding to a single free phase 
among the output fields. 

To solve the pump equations we first extract the 
dependence of the phases of the signal and idler on the 
phases of the pump, as discussed in section IIII Bl This 
gives 

e iargp t = |^| e iargF T) (10) 

where Lf is the left-hand side of Eq. [3] evaluated for 
argp| = 0. Taking the modulus of Eq. E3 w e have the 
general form 

\a+\z\ 2 f3\ = \F,\, (11) 

where a and (3 are functions of the intensities of the pump 
fields. We solve Eq. ^]to determine the output intensi- 
ties that are consistent with the strength of the pump-up 
driving, as functions of the intensities of the pump fields 
if (if ) and if . We then solve the spin-flipped version of 
Eq. llll to determine the output intensities consistent with 
the strength of the pump-down driving. Equating these 
two intensities gives a nonlinear equation which we solve 
to determine if as a function of the external pumps. 

Fig0] illustrates the solution to the pump equations 
(0. The damping and Hopficld coefficients for the sig- 
nal and idler are as in Fig^ For the pump we have 
used 7 P = O.lmeV and \x p \ =0.8, which we again es- 
timate to be appropriate to the system reported in Ref. 
HH We have taken an external drive of fixed circular- 
ity, a = 0.1, and varying intensity. The top panel shows 
the intensities of the pump fields, and the bottom panel 
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shows the output intensity. Increasing the driving inten- 
sity from zero we first find two continuous thresholds, 
where steady-states appear starting with zero output in- 
tensity. The pump fields at these thresholds are marked 
as filled dots on Figd At these points, the pump inten- 
sities match, up to a common factor, with the driving 
intensities: /jj = |i*\jj 2 /| A p | 2 . For this small value of 
g, they are approximately the thresholds for each polari- 
sation of the uncoupled device. With increasing driving, 
the output intensity in each of these steady-states in- 
creases and, in constrast to the uncoupled device, the 
pump fields change. Increasing the driving intensity still 
further, we find a third threshold at which a new steady- 
state appears, and then splits into two states. This third 
threshold is discontinuous, i.e. the solution appears with 
a finite output intensity. It corresponds to the pump 
fields marked with the open dot on Fig. ^ 



IV. DISCUSSION 

In Ref. UJJ, the intensity and circularity of the signal 
output was measured under continuous driving of the 
pump field. The driving circularity was varied from cir- 
cular to linear and the total pump intensity was of the 
same order of magnitude as the threshold for a circularly 
polarised pump. As the pump circularity a was decreased 
from one, the signal intensity increased by a factor of five 
to a maximum at a ss 0.4, and then decreased again. The 
circularity of the signal was approximately constant, ex- 
cept near to linear pumping, where it dropped to zero. 

Considering the pair process underlying the paramet- 
ric oscillator, one might expect that without polarisa- 
tion coupling the output would monotonically reduce 
as the pump circularity is reduced: scattering only oc- 
curs within each polarisation, and moving away from cir- 
cularity reduces the population of each polarisationAi. 
However, such an interpretation overlooks the effects of 
pump depletion and coherence. Because of these ef- 
fects, the steady-state output of a single polarisation 
pumped with intensity i+ is actually proportional to 
\fl+ — V ^thresh (Ref. 0). For total driving intensities 
greater than twice the single-polarisation threshold, there 
is a critical pump circularity below which both polarisa- 
tions are above threshold. Below this critical circularity, 
the total output increases as the pump circularity is re- 
duced, with a local maximum for a linear pump. Thus an 
increase in the output as the pump circularity is reduced 
does not in general imply the existence of polarisation 
coupling. However, without such a coupling it is difficult 
to explain the de pen dence of output on pump circularity 
reported in Ref. llOt the enhancement away from circu- 
larity is too strong, and the output peaks for elliptical, 
rather than circular or linear, pumping. 

The lower branches shown on Fig[5]have no structure 
to suggest that they would give a maximum output for an 
clliptically polarised pump. However, the upper branches 
illustrated for g = 5 and g — 7 do have such structure. 



As the pump circularity is reduced from one, we would 
first go up through the threshold for these solutions, then 
move back towards it, and for some parameters go below 
it again as we approach linear pumping. The circularity 
of the turning points illustrated for g = 7 is a = 0.47. 
While this is roughly consistent with the experimental re- 
sults, a detailed fit to the experiment is beyond the scope 
of the present paper. It would involve a large number of 
parameters, and we expect the results to be sensitive to 
details left out of the present model such as the blueshifts. 

It is only stable steady-states which are relevant to 
continuously pumped experiments. We have analysed 
the stability of some of the steady-state solutions to our 
model for a special case in which all the effective damping 
rates are equal, 7^ = j' s = 7-. For g = 0.1, there are 
two solutions with continuous thresholds and two with 
discontinuous thresholds, as there arc in Fig.^J We find 
that only the solution with the lowest threshold is stable. 
For g = 2.0, we find only one solution with a continuous 
threshold, as well as two with discontinuous thresholds. 
In that case, both the continuous solution and one of the 
discontinuous solutions is stable. Thus it is possible for 
our model to have stable solutions other than that with 
the lowest threshold, and to have more than one stable 
solution. 

The pulsed experiments of Ref. ^] have recently been 
addressed by Kavokin et Their theory reproduces 

the experimentally observed rotations of the linear polar- 
isation without scattering between polarisations. In their 
theory, the rotation of the linear polarisation comes from 
the different blueshifts of the two circular polarisation 
states. The present work does not include this effect, be- 
cause we have taken the detunings of the two polarisation 
states to be the same. We expect that if the blueshifts 
are small compared with the interactions between polar- 
isations then the steady-states will only contain a single 
frequency for the pump, signal and idler, and the treat- 
ment given here will be qualitatively correct. However, a 
splitting of the two circular polarisation states might be 
an alternative explanation for the steady-states results of 
Ref. US 

The agreement between the polarisation rotations seen 
in the pulsed experiments^ and the theory of Ref. sug- 
gests that a polarisation coupling term is not required 
to explain these results. However, this does not imply 
that the polarisation coupling is always irrelevant. The 
polarisation rotations produced by a coupling could be 
small compared with those produced by the blueshifts, 
allowing a good fit to this aspect of the data without 
a coupling. Note also that the pulsed experiments are 
done at much higher excitation powers, typically around 
a hundred times greater, than the steady-state experi- 
ments. The polarisation coupling could also depend on 
sample parameters such as the energy difference between 
the pump polaritons and the biexcitonifi. 

An interaction of the form we propose corresponds mi- 
croscopically to an interaction between excitons of oppo- 
site spin, which could come from the Coulomb intcrac- 
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tions between the electrons and the holes in the excitons. 
To first order in these Coulomb interactions, the inter- 
action between opposite spin excitons is negligibl o 1 ^ 1 ^ , 
because in that approximation small wavevector scatter- 
ing is dominated by electron-electron and hole-hole ex- 
change. Thus a finite value of g implies the significance 
of higher-order processes, which can produce interactions 
between opposite spin excitonsi^. We suggest that the 
higher-order processes could involve the rrij = ±2 exci- 
tons, which are close in energy, and produce an interac- 
tion at second order. Another recent suggestion^ is that 
the interaction involves excited states of the excitons. 



V. CONCLUSIONS 

We have studied the steady-states of a model of micro- 
cavity parametric oscillation with coupled polarisations. 
For small values of the coupling, we find two steady-states 
corresponding to those of the uncoupled device. How- 
ever, due to the increased scope for arranging the pump 
depletion, there are also two steady-states which appear 
discontinuously, i.e. with a finite value of the output in- 
tensity, as the driving intensity is increased. For general 
values of the coupling Fig|3 suggests that there will be 
cither one or two continuous solutions, depending on the 



coupling and pump circularity. There may also be dis- 
continuous solutions. For some parameters more than 
one steady-state can be stable, in which case it should be 
possible to observe switching between the states induced 
either by noise or by external probes. 

The coupling between polarisations introduces two 
types of mixing term into the equations of motion for 
the fields. As well as the straightforward analog of the 
process considered in the uncoupled model, there are pro- 
cesses which exchange the spins of two fields. Such spin- 
flip processes lead to a rotation of the output polarisation 
with respect to the pump. They also produce shifts in 
the output frequencies with varying pump intensity or 
circularity, even in the absence of the shifts associated 
with the mean-field exciton-exciton interaction. 
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